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ABSTRACT 

There is much disagreement as to (a) what causes clear air 
turbulence (turbulence which is not in or near convective clouds and 
is above 15,000 feet in altitude) and (b) which Wolcer olepiical param- 
eters can be used to detect and forecast its occurrence. The approach 
to this problem has been to relate not one parameter to clear air 
turbulence but various parameters. By summing eae parameters areas 
can be defined where there is a high probability of encountering clear 
air turbulence. Each paramete fHGSDORnR s ed on a statistical study 
which found a relationship with clear air turbulence. The parameters 
used were horizontal and vertical shear, curvature, kinetic energy 
and their derivatives. 

The numerical forecasting program proposed here can be extended 
to the stratosphere when more reliable height and temperature fields 
are available. This program will have much more significance when 
intermediate forecast height fields, temperature fields and a grid of 


much smaller mesh length are available. 
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i. INTRODUCTION 

The phenomenon of clear air turbulence (henceforth denoted as CAT) 
appears to consist of random three dimensional eddies that occur in 
certain confined regions of the atmosphere. This phenomenon was first 
recognized in the early 1940's as “air pockets", and has gained in im- 
portance and depth of study with the development of fast-flying swept 
wing and delta winged aircraft. It is of prime importance to the avia- 
tion industry because it affects the safety and comfort of passengers 
and crew, as well as operational costs. There have been several cases 
where loss of control of aircraft, structural damage, passenger and 
crew injuries (even deaths) have resulted from CAT. 

Aircraft manufacturers and the airlines are interested in CAT in 
order to determine the strength of airframe design so that it can be 
built to withstand all degrees of turbulence. Planning of supersonic 
transport aircraft is considering the effects of this unexpected turb- 
ulence occurring at any flight level in the atmosphere. Therefore, it 
becomes an even greater operational problem in the era of future design 
of SST aircraft as to cost, safety and comfort. CAT is usually less 
intense than turbulence encountered in thunderstorms 1 . CAT may be 
considered at times more dangerous than "thunderstorm turbulence” as 
it usually occurs with no visible warning. 

The various military agencies are deeply concerned with the under- 
standing of CAT and its prediction so they may be better prepared to 
take it into account operationally in all areas of the world. Scientists 
are highly interested in CAT because it is a phenomenon of our environ- 
ment that is not clearly understood. We are faced with the problem of 
having to rely on mesoscale data while the phenomenon is of a microscale. 


Because of the great differences in scales involved in the forecasting 


problem, isolated patches of CAT (in either space or time) are very 
difficult to identify. It is the intent of this paper to develop a 
numerical forecasting system of CAT whereby an empirical approach has 


been used to accomplish this end. 


2. DATA COLLECTION 

The nature of clear air turbulence, its physics and its meteorology, 
are still not completely understood. More mesoscale (2km to 100km 
horizontal distance) and microscale (less than 2km) studies must be 
made concerning the causes and generation of CAT. Any improvement in 
forecasting skill that may be realized by this study will probably result 
from increased mesoscale input data into the forecast problem. The only 
mesoscale measurements available operationally at this time which concern 
the problem are rawinsonde data that give a nearly continuous measure- 
ment in the vertical of wind, pressure, temperature and humidity from 
the surface to above 30km in altitude, and satellite cloud pictures 
that delineate cloud structures down to a scale of about 3kn. 

The usual definition of CAT is "atmospheric turbulence which is 
not in or near convective clouds, including thunderstorms, and is not 
below 15,000 feet in altitude"'. (4) Thus, mechanical turbulence induced 
by rough terrain is excluded. It is realized that this is an extremely 
arbitrary definition resulting from a desire to simplify pilot reporting 
procedures. 

Turbulence intensities are, at present, designated as light, mod- 
erate, severe and extreme. In recent years several definitions of the 
four categories of turbulence have been proposed. The current defin- 
itions in official use were developed by the National Advisory Committee 
for Aeronautics (NACA) Subcommittee on Meteorological Problems (1957). 


See Table l. 


There are few quantitative measurements of atmospheric turbulence 
at any scale. The bulk of the aircraft turbulence data available is 
based upon the four categories stated in Table 1, (4) ; and, conse- 
quently is highly subjective and qualitative in nature. Included in the 
factors that affect the subjective decision of the pilot when reporting 
turbulence are the wing loading, the aircraft's speed and aes the 
pilot's training, experience, and mee psychological reactions (1) . 

Objective criteria for turbulence reporting must wait until more 
measurements are made of the conditions existing in the turbulent eddies, 
and until cockpit instrumentation includes a gust load or turbulence 
indicator. 

Much time, effort and financial support have been expended in this 
country by the government, the aircraft industry, and private and public 
research facilities to investigate the problem of clear air turbulence 
and to find a good approach leading to a solution to this problem. This 
research has gone forward in three broad categories: first, research 
concerning the correlation of CAT with macroscale, mesoscale, and/or 
microscale atmospheric measurements; second, research into instrumen- 
tation for the detection of CAT sufficiently far in advance of the air- 
craft to allow evasive measures (6) and third, research by military and 
airline organizations concerning the operational aspects of CAT. 

Various projects such as HICAT, ALLCAT, and TOPCAT have been under- 
taken to study clear air turbulence. The problems encountered were many 
but the results showed that it was indeed possible, with suitable in- 
strumentation, to find, track and record CAT. The use of constant level 
ballons and doppler radar to detect CAT has met with some success, but 
again the lack of detail necessary in CAT studies leaves much to be 


desired. 
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Measurements of the microstructure, which contains the perturbations 
of CAT dimensions (100 to 500 meter wavelength), have been taken by 
specially instrumented aircraft. Until recently the aircraft itself 
was used as a sensor to measure atmospheric gusts from the aircraft 
acceleration data. The Air Force realized the danger of relying solely 
on the acceleration response of the U-2 aircraft as a measure of turbu- 
lence at high altitudes. A full knowledge of the aircraft's response 
to turbulence over a wide range of wavelengths is required for meaning- 
ful interpretation of such data. This procedure usually gives reliable 
results at short wavelengths up to a few hundred meters. At long wave- 
lengths, this instrumentation becomes less sensitive. Accelerations in 
longer waves are usually small and can be masked by pilot induced air- 
craft motions. 

A supersonic or hypersonic craft of some radical shape, flying four 
to ten times the speed of the U-2, will obviously have a somewhat dif- 
ferent response to turbulence than the U-2. An aircraft flying at 
these high speeds would be affected much more by longer turbulence wave- 
lengths and less by the shorter. 

As pointed out earlier, the meteorologist has few direct measure- 
ments of turbulence intensities and must depend upon the accuracy of 
the intensities reported by pilots. In our study, we used the Colson 
monthly turbulence summaries which were obtained from the Air Force 
3rd Weather Wing at Offutt AFB, Omaha, Neb. This report suited our 
needs most closely because the CAT reports were detailed as to location, 
time, altitude and intensity. In the period used in our study, December 


1964 through March 1965, there were no less than 3670 CAT reports. The 
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reports were located over the United States and extended from about 
15,000 feet to over 45,000 feet. The data were collected from military, 
civilian and private aircraft. They are, as mentioned before, quite 
subjective. 

In the course of this paper three time periods or intervals will 
be used. It is important to establish at this point which periods were 
used and for what purpose. 

A four month period (December 1964 through March 1965) is used for 
the research program during which several fields were constructed. 
Patterns of fields of different meteorological parameters were some- 
times similar. Then all but one of them were dropped in our further 
investigations. In this way only three of the original seven parameters 
were retained. 

The three day period (10 - 12 March, 1965) was used for a verifi- 
cation or correlation study. During this period of time a comparison 
was made to find out how many times CAT actually occurred in areas where 
it was predicted. 

Finally, 23 February 1965 was arbitrarily chosen from the four 
month period and used only for illustrating the various fields and 


parameters used in the research program. 
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3. THE CLEAR AIR TURBULENCE STUDY 

Clear air turbulence (CAT) is a microscale phenomenon (less than 
2 kilometers in horizontal extent) but the conditions which are sympto- 
matic of its existence are of snyoptic scale. Therefore, synoptic 
scale parameters may be used to determine areas where CAT could occur, 
that is forecasting areas where there is a high probability of CAT. 

From various reports and their contradicaions it appears that no 
single parameter can detect CAT. Certain parameters can detect the 
possibility of CAT in some synoptic situations but fail in others. 

The approach used in this paper was to take statistical studies 
made by previous investigators and to use the basic parameters which 
they related to CAT (2,3,5,7,8). If the magnitude of any one parameter 
becomes large or the sum of various parameters becomes large then there 
is a high probability of CAT in that area. Therefore, the problem is 
not one of forecasting actual CAT but rather to forecast areas of high 
and low probability of encountering CAT. In this way flights may be 
planned so as to expect least CAT. 

The research program consisted of computing the equations shown 
in Appendix A. These equations were computed using the CDC 1604 digital 
computer. The program was written utilizing symbolic coded relocatable 
assembly program (SCRAP). It was necessary to use fixed point fractional 
numbers in order that Fleet Numerical Weather Facility (FNWF) subroutines 
could be used. All finite differences were computed using standard FNWF 
mesh length of 381 kilometers true at 60 degrees north latitude. There 
are no time derivatives in the program. The research reported here is 


accomplished using analytical fields. Operational use would employ 
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forecast fields. Resulics would natvraily noc hav- peer ase ~ 
cast fields wermlused in this reswarch. ‘Prog 24 hous’ printydé 2f whe 
bottom of each tield represents the practic2l forecast wnter alco oe 
used operationally. 

The research program was written to compute CAT in three layers 
500 to 300, 300 to 200, and 200 to 100 millibars. Computations were 
not made for the third layer since 100 millibar fields were not available. 

The research program was written to compute on the entire 63 x 63 
FNWF grid of the northern hemisphere. A boundary condition of zero was 
usea for the outside rows and columns. The print routines are 22 x 22 
extracts of the United States starting at FNWF grid point J008, I018. 
The latitude and longitude coordinates of the four corners of the printed 
fields starting with the lower left corner proceeding clockwise are: 

9 As OO lh SN Gol 7 IN. 3 alleles ON BS ll, 

All printed fiele- are pure numbers and have no dimensional 
meaning. All scaled outputs have been shifted to the left end of the 
register and the fFl-st three numbers with sign bit are printed out in 
decimal. Crid points can take on values from -999 to +999 except those 
fields which have been made all positive. The decimal point does not 
appiar on th» vcrinted fields. Therefore, the printed grid point values 
are from -999 to +999, 

In the research program the capital letters refer to the parameter 
as computed from 3 data. The small case letters serving as exponents 
represent the number of times and direction the register has to be 
shifted in order to place the significant portion into the first three 
numbers. Therefore, the two |vith exponents represent the scaling co- 


efficient. Since the cemputations were in fixee hoint *-actione! al) 
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printouts had to be shifted so as not to exceed one at any grid point 
in the field. Exceeding one would result in a meaningless value at 
that particular grid point. In addition the entire field had to be 
kept large enough so that patterns could exist and not be at or near 
zero throughout the field. Since these fields are summed they must be 
small enough so as not to cause the summation field to exceed one at 
any grid point. 

The research program was run for thirty-six days during December 
1964 and January, February, and March 1965. The thirty-six days were 
chosen because they were the most active in CAT reports during the four 
month period. In other words, there were more reports by pilots en- 
countering CAT on these particular days. In order to show an example 
of each field printed out by the research program the time 00Z 23 FEB 65 
was arbitrarily selected. These fields appear in Appendix C. Each field 


was produced by an individual term which will be described as follows: 


THE FIRST TERM OF THE RESEARCH PROGRAM (APPENDIX A) 

The First Term is 2° A where A is the Laplacian of absolute 
vortisity. When this term is negative there is a local maximum of 
absolute vorticity meaning it is larger at that grid point than the 
average of the surrounding grid points. Therefore, the cyclonic curva- 
ture or cyclonic shear or both are relatively large at that grid point. 
This should correspond to the cold side of the jet especially in 
throughs. According to Endlich and McLean (1) there is a greater per- 
centage of CAT on the cold side of the jet. Also according to Harrison 
(2) there is a strong tendency for moderate to severe CAT to be asso- 


ciated with trough lines. 


BS 


This term was computed for the layer by first calculating the 
absolute vorticity of the upper and lower level D fields. The Laplacian 
was then taken of each field and a vertical average made of the upper 
and lower levels to obtain the Laplacian of vorticity of the layer. 

This parameter has depicted most of the CAT associated with the 
trough over the western United States. The severe CAT near New Orleans 
is in an area of very large negative numbers. The field is contoured 
at intervals of 100 with the origin at zero. 

THE SECOND TERM 

The Second Term is 2? B where B is the absolute value of the 
vertical change in the vector thermal wind. One of the parameters 
which Lake's (7) statistical testing indicated was associated with CAT 
was the vertical gradient of wind shear. As shown by Richardson (9) 
the thermal wind shear is proportional to the gradient of static 
stability. 

The u and v components of the thermal wind were computed from the 
upper and lower level temperature fields. The difference between the 
upper and lower level values of the u component was found and each 
difference was squared. This was also done for the v component. The 
Square root of the sum of the squared differences gives the magnitude 
of the vector difference. According to Endlich and McLean (3) the 
largest values of the thermal wind shear appear on the warm side of the 
jet. This was found to be true throughout the four months. The contour 
interval for this field is 25 and the origin is zero. 

THE THIRD TERM 

The Third Term is 2° C where C is one half the geostrophic wind 

velocity squared. C therefore represents the specific kinetic energy 


or in other words the kinetic energy per unit mass. 
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Clem (2) found that most cases of moderate to severe CAT were 
associated with areas of isotach maxima. 

This term was computed for the layer by calculating the u and v 
components of the geostrophic wind at the upper and lower levels. The 
upper and lower level u components were vertically averaged to obtain 
an average u component for the layer. The average v component for the 
layer was obtained by a similar process. The magnitude of the velocity 
squared is just the sum of the squared components. 

In the research program this field is contoured at intervals of 
100 with the origin at zero. Contoured at this interval the kinetic 
energy field clearly depicts the isotach maxima regions. The kinetic 
energy field in Appendix C shows this field depicting the CAT in the 
western part of the United States occurring in areas of relatively 
large wind speeds. The kinetic energy field fails to indicate the 
severe CAT near New Orleans because it occurs in an area of relatively 


light winds. 


THE FOURTH TERM 

The Fourth Term is 2¢ D where D is the absolute value of the 
derivative of the kinetic energy with respect to pressure. 

Lake's (7) statistical testing indicated that the gust intensities 
are related to the vertical gradients of horizontal kinetic energy. 
This term was computed for the layer by first calculating the velocity 
squared at the upper and Lower levels. The vertical gradient for the 
layer was then obtained by computing the difference between the upper 
and lower level values of the velocity squared. The absolute value was 


taken so as to have all values positive. CAT should be associated with 
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large values of this field. Large values of this field were found 
only in areas of large values of kinetic energy. This field was 
therefore redundant and was eliminated from the CAT forecast program. 


The contour interval was 100 with the origin at zero. 


THE FIFTH TERM 

The Fifth Term is 2° E where E is the absolute value of the 
Laplacian of kinetic energy. The statistical survey made by Endlich 
and McLean (3) shows the maximum occurrence of CAT along the edges of 
the isotach maxima. The Laplacian of kinetic energy shows large 
horizontal changes in kinetic energy, both positive and negative. 
Therefore the absolute value of the term is taken in order to give 
only positive numbers. The contour interval was 25 with the origin 
@ee ae LO, 

This term was introduced to depict the areas of large horizontal 
change in kinetic energy. However, there was no relationship with 
reported CAT. This field was therefore eliminated from the CAT fore- 


cast program. 


THE SIXTH TERM 

At this point in the research program it was necessary to sum the 
first five terms due to computer memory space. This term, referred to 
as KAT1, was the Sixth Term. The contour interval was 250 with the 
origin at zero. It was found that this term did not supply signif- 
icantly new information since it was dominated by the kinetic energy 


and the two associated terms. 


THE SEVENTH TERM 
The Seventh Term is 2' F where F is the Jacobian of temperature 


and omega (component of the wind normal to the pressure surface). This 
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term was developed by Dr. Moore of Douglas Aircraft and Dr. Krishnamurti 
(8). The latter was associated at that time with the University of 
California, Los Angeles and consultant to Douglas Aircraft. The term 
was developed as the Jacobian of temperature and three dimensional 
divergence. As shown in their paper this is proportional to the nega- 
tive of the Jacobian of temperature and omega. This term was computed 
for the lower level of each layer in the research program. The contour 
interval for this term was 100 with the origin at zero. In this pro- 
gram no significant relationship was found with large negative or 
positive numbers over the four month period. This term was therefore 


eliminated from the CAT forecast program. 


THE EIGHTH TERM 

The Eighth Term is 2&8 G where G is the absolute value of horizontal 
divergence. It was computed by taking the derivative of omega with 
respect to pressure. As previously stated the entrance and exit regions 
of isotach maxima areas have been found to be associated with CAT. 

These areas are also associated with horizontal convergence at the 
entrance and divergence at the exit regions. Therefore the absolute 
value of the change in omega with pressure represents the convergence 
and divergence in the layer parallel to the pressure surfaces. 

This term was computed by subtracting the lower level omega value 
from the upper level omega value at each grid point. Areas of con- 
vergence and divergence of the height field are quite vividly depicted 
by the divergence field. No significant relationship was found between 
the divergence field and the CAT occurrences, therefore, it was elimin- 


ated from the CAT forecast program. 
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THE NINTH TERM 

The Ninth Term was KAT2,the summation of all previous terms. No 
significant relationship was found with CAT occurrences since several 
terms tended to cancel each other out. 

The program then goes into the second layer from 300 to 200 milli- 
bars. All terms were computed and scaled the same, except the diver- 
gence term. It could not be computed because the 200 millibar omega 
field was not available. 

The third layer from 200 to 100 millibars could not be run for 
these time periods since the 100 millibar fields were not available. 
There were very few CAT reports above 200 millibars, therefore, the 
loss was insignificant. 

In all three layers the lower level height field is printed out 
first. This gives a general impression of the synoptic situation and 
renders more significance to the patterns developed in the other fields. 
The contour interval for the 500 millibar field is 60 meters with the 
origin at 5580 meters. The contour interval for the 300 millibar field 
is 120 meters with the origin at 9120 meters. The contour interval for 
the 200 millibar field is 120 meters with the origin at 11,760 meters. 

The CAT forecast program appears in Appendix D. The first three 
terms of the research program are used with minor changes. In the first 
term "a'' is changed to minus one and the contour interval has been 
changed to 150 to give better defined patterms. The second term has 
been used unchanged. The third term is unchanged except for the con- 
tour interval which was changed to 50 to increase the pattern size. 

The KAT field itself is the summation of these three terms and gives 


quite reasonable patterns and pattern sizes. 
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The pattern area depicting a high probability of CAT would neces- 
sarily be larger during a more active CAT period. The most active 
part of the year was the four month period December 1964, January, 
February, and March 1965. During these months the most active three 
day period was the tenth through the twelfth of March 1965. Therefore 
the KAT fields have quite large pattern sizes in Appendix E since they 
represent the most active three days of the year. The KAT field 
patterns were smaller for less active periods. The total area covered 
by these patterns is much less in the KAT field than in the other three 
parameter fields. 

This is exactly what was attempted in order to obtain optimum 
size of the forecasted danger areas. If the KAT field patterns are 
too large, flights will be rerouted unnecessarily. On the other hand, 
if the KAT field patterns are too small, there is a real danger of CAT 
occurring outside these areas. Therefore, the restraint of the KAT 
field patterns is necessary in order to have an operationally useable 


product. 
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4, DISCUSSION AND RECOMMENDATION FOR FUTURE STUDIES 

The period used in this paper was chosen because of the largest 
number of reported CAT occurrences. Of the four month period (December 
1964, January, February and March 1965) there were scattered periods 
where a large number of occurrences were reported. This four-month 
period was used to determine which parameters were best suited for 
forecasting purposes. The three-day period of 10-12 March 1965 was 
selected for a correlation study in order to find out how successful 
our forecast method is. Tables 2, 3, and 4 show the various fields 
used and the resultant KAT field for the 10th, 11th and 12th of March 
1965. 

The use of the term "percent correlation'' as used in this study 
does not mean to imply a statistical correlation. The ideal forecast 
verification makes use of those cases where CAT is forecast, but does 
not occur, and where CAT is not forecast and does not occur. In our 
study it was impossible to take those cases quantitatively into account. 
Therefore, it is to be understood that ''correlation'’ as used in this 
study was a general comparison of those reported CAT occurrences that 
fell within the delineated area of high probability of CAT against those 
that did not. For example, if there were ten reported CAT occurrences 
for a given layer and time period and six of these reports fell inside 
or on the line delineating the CAT area and four reports fell outside 
the area, then for that field, layer and time period we would List it 
as six occurrences correlated or a sixty percent correlation. 

Listed are the names of the fields, the number of occurrences 
of CAT for each field and the percent correlation by field. Also shown 
are the number of occurrences and the number of occurrences that cor- 


related by field and CAT intensity. 
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After combining the Laplacian of Vorticity, Vertical Gradient of 
Thermal Wind and Kinetic Energy we arrive at the KAT Field which is our 
end product for the area of high probability of CAT occurrences. Even 
though a higher correlation may be seen in some fields other than the 
KAT field, one should realize that these fields encompassed a larger 
than average area. In such cases one must expect a high correlation. 

The high correlation is then not due to the finesse of the forecast 
method, but rather due to the fact that for most of the USA there was 
a forecast of high probability for CAT. Theoretically, it would be a 
good idea to divide the percent correlation (such as we computed) by 
the size of the area for which CAT was forecasted. We did not follow 
this idea quantitatively, but only qualitatively. Therefore, one finds 
that the percent correlation for our ultimate forecast (labeled KAT) 
is sometimes lower than the percent correlation for one of the three 
separate forecasting fields. 

Table 5 is a summary for the three day period. It shows the total 
number of CAT occurrences by turbulence category, percent correlation 
by field and turbulence category, and the three day percent correlation 
by field. 

The results were most encouraging and we feel that our end product 
was a substantial step in at least the right direction toward forecasting 
clear air turbulence. Our knowledge of the meso- and micro- structure 
of flow patterns in the free atmosphere, especially above the tropopause, 
is still rather poor. A strong need still exists for a well organized 
and well equipped measurement program, especially at flight levels of 


the future supersonic transport aircraft. Measurement programs using 
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methods of data collection other than aircraft should be sought in order 
to obtain more information on the real micro- structure of the atmosphere, 
without the large disturbances which a flying aircraft will create 

itself. 

Case studies of CAT occurrences so far were limited to a comparison 
of turbulence location with atmospheric parameters measured as closely 
as possible to the time of occurrence. We might gain some additional 
information on the physical causes of CAT if the development and pre- 
vious history of flow patterns bearing CAT were studied. 

There still exists a need for sensitive, accurate and compact 
instrumentation, especially an accelerometer which measures and records 
the three components of gustiness separately and simultaneously. 

In summarizing we would like to state that turbulence research in 
the free atmosphere has come a long way, especially when we consider 
the fact that measurements are most difficult to duplicate in controlled 
experiments. We need the free atmosphere above us to conduct our re- 
search, and this same atmosphere has an infinite choice of parameter 
combinations. This should provide the seed of interest to combine the 
efforts of physics, aerodynamics, mathematics, statistics and meteor- 
ology to seek out and find a more complete and more satisfying solution 


of the problem. 
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The numbers in parentheses on this and the following pages refer 
to the page numbers in the text where each item is discussed. 
The numbers below each term in the following two equations identify 


the number of each term (e.g., KAT1 is the sixth term). 
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9. APPENDIX D 


THE CLEAR AIR TURBULENCE FORECAST COMPUIER PROGRAM 
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TABLE 1 


AIRCRAFT TURBULENCE CRITERIA 


CATEGORY DEFINITION 
LIGHT A turbulent condition during which 


occupants may be required to use 
seat belts, but objects in the 
aircraft remain at rest. 

MODE RATE A turbulent condition in which 
occupants require seat belts and 
occasionally are thrown against the 
belt. Unsecured objects in aircraft 
move about. 

SEVERE A turbulent condition in which the 
aircraft momentarily may be out of 
control. Occupants are thrown 
violently against the belt and back 
into the seat. Objects not secured 
in the aircraft are tossed about. 

EXTREME A rarely encountered turbulent con- 
dition in which the aircraft is 
violently tossed about, and is 
practically impossible to control. 
May cause structural damage to the 


arceract. 





163 


¢ Aldvai 


es 
ae ROWING DILANTY 
06 , en GNIM TVWYSHL JO INFIGVYD TWO1LYGA 
6 CT a eee ALIOILYOA JO NVIOVTAV1 


vet Ee il 
1 ATT 


re 
TET am 


Led GNIM ‘IVWYSHL JO INSLIVYD IVOILYZ, 


CQO YWN Ol ALIDILYOA JO NVIOVId‘1 
+. Ye eo qidia .vy 


cay 
mH 
uw) 


0 


ON 





wy) wy) 
— rd 


L 








ALIOILYOA JO NVIOVITdVI 


wy) 
ur) on o * or 
= ~— 


N N 
N N 
sO \O =. fm ot 
OG) ™ f it Oo 


ON} rm 
ae 


= 





N 
™I] 
Cl 
ON 








OYaNd OILANIN | 





Cee 
OE 
ea aad rieg | een 


oo HOD FAOOO a a YdOOFAOOO] YYOD FAYDOO AO 
op - 5 7: df: 











4: TV LOL 





164 


f / 





ce 
aoe 
Pe 
re 
ra 


ep et ' 
41 
co 
—t 
Sa 


rt" 


Pe 
Pr 
oa 
os 
ahi 








i 


re 


ere 


ddOo 


conn 


~t 
— 
© S © 
Ww \O cO 


on 
~—a 
ra 
ts 
— 
f™ 
tf 
—! N 
N wy 


~t 
— 


ete 
eT 
ca 
BoE 
Ts 
Pies 
GE 
rhe 


€ ATavd 


ca 
t 


~ 


ate 
Oo 


O 
CO 


Ww) 


© 
as 


— Re nl ke 
bE EI 
a 


8 
O 
& 


SAA ‘T 


- 
2 
= 
© 
= 


és 





S 


AW LL 


G9 UV IT 


Ww 


9 YVW [I 


G9 UVW IT 


CO UVN Wt 


Lit 


Uaala LV 


WOYXANaA O1 LAINIA 


CNM video menerLivdo TV LAaA 


iota A NVIOVId¥ 1 


Cat tev 


medava JLLINI a 
GQNIM TVWNYSHL JO INAICVeS TVWOT LSA 


ALIOILYOA JO NVIOVIdV' 
Cererre ay a 
AQYANA OILANTIO 


GNIM ‘IVWHSHL JO INAICVaD ‘TVOILYdA 


ALIOTLAOA dO NVIOVIdV' 


QTAT deel V2 


AOYANA DILANIN 





CNIM TVWWHAHL dO LNAICVeS ‘TVOLLddA 


ALIOILYOA JO NVIOVIdV'l 


ed yal 











la 


7 ATAVL 


ovata LVI 


AOQYHNd OL LININ 


QNIM ‘IVWHHYHL AO INSIGVYS TWOILYHA 


rei 
— 
tH 
NY 
GN 
em 


ALIOILYOA JO NVIOVTdVI 


rita Lv 


AQYYNA OI LANIN 


x se) rand 
aH = oo 
CO N — ~t ~t ees 
CO ON CO ON ON CO 


\O Oo 
N 





—t 
re 


© 
7 
_—i 


QNIM ‘TVWHSHL dO LNSIGVYS TVOILYAA 


wy) 
—| 


Zo So ov cl ALIOILAYOA AO NVIOVTdV'! 


wtae LV 


AOYANA OILYNIA 


— 








| 1 = a NS N N NS 
-_ —s - — 
= ~) 
—_— 
st mals ~t wy un ny) Fa) N N N N 
—=| asa ei — i 4 = — 


GQNIM ‘IVWHAHL AO LNAITGVYS TWOILYHA 


fede] ef bel et of of op - 
er — — 


— 


eH 
b— 


61 ZOO 7S9 UYVW Cl Ailowrden AO NVTOV dV 


On CO; rm 
\O +] ME 


tat 
ee — 
ee ed tet re SE 
Ttholn tenth 1 te LL Ak. (NIM TVWYSHL JO LINAIGVYD TVDILYSA 
Riese ieitedeieteit It) ee 


yYuoo] WOOT uNOD Tt. YOO VYODTYODO!] yyODTyDD0 AO 
4: 4: d: 4 4: d: J: 7 YAAVTE ANIL qLvd SQTdaI4 


uw) 
\O 
~% 
= 
N 
rj 











166 


GC WTAV.L 


eo 
Oe 
O06 


COMM AT wn ddH LO INITIVYD AVORIYAA 








ALIOILYOA dO NVIOV1dVI 





( GTAI4 Aq ) 





NOLLVIAYdOD LNAOYAd 








O9 


c 


oy = 
C7) try 
oO 
™ 

—! 

mM ™ ™ 


NS 
™ 


fe 


© 


0 


a 
Ya) 
ON 


GNIM ‘IVNYAHL JO LNAIGVYS ‘TIVOILYAA 


\O 
= 
™ 


oO 
a 
tf 


C9 ALIOILYOA JO NVIDVIdVTI 


( KYOODALVO GNV GIAIA Ad ) 


NOMmLVidddO) INIDudd 





(C961 HOYVW ZI OL S96T HOV OT) 








oO 
NS 
~) 


0 


om 
uy 





om COE SHONAYANDDO IV JO YaAGWNN 'TIVLOL 


aaa AAS 
Ol 


ALVaAGOK 


ha GA aS ALVaa don 


OL SatTi08aqe) 


ADVYed don LHOTT LIV 











105 


Lie 


123 


INITIAL DISTRIBUTION LIST 


No. 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library 
Naval Postgraduate School 
Monterey, California 93940 


Professor William van der Bijl 

Department of Meteorology and Oceanography 
Naval Postgraduate School 

Monterey, California 93940 


Office of the Naval Weather Service 
Naval Station (Washington Navy Yard Annex) 
Washington, D. C. 20390 


Department of Meteorology and Oceanography 
Naval Postgraduate School 
Monterey, California 93940 


LT William A. Morgan 
4500 Curtiss Drive 
Virginia Beach, Va. 23455 


LCDR Michael J. Ettel 
*% VF 101 

NAS Oceana 

Virginia Beach, Va. 


Mr. J. Ehernberger 

NASA Flight Research Center 
POR Boxe 273 

Edwards AFB, California 93523 


Mr. Walter M. Crooks 
Lockheed California Company 
Burbank, California 


Dr. R. L. Moore 

Douglas Aircraft Company 
Aircraft Division 

Long Beach, California 


CAPT De Giovanni 
3rd Weather Wing (3V) 
Offutt AFB, Nebraska 68113 


Physics Department 
ot. Johns University 
Collegeville, Minnesota 


Copies 


20 





Jey 


14. 


by. 


ro’. 


7s 


lier. 


ie 


20. 


Dike 


Zee 


Z3- 


24. 


No. 


Professor R. J. Renard 
Naval Postgraduate School 
Monterey, California 93940 


Professor T. N. Krishnamurti 
Naval Postgraduate School 
Monterey, California 93940 


Officer in Charge 

Naval Weather Research Facility 
Naval Air Station, Building R-48 
Norfolk, Virginia 23511 


Commanding Officer 
Fleet Weather Central 
Navy Department 
Washington D. C. 20390 


Officer in Charge 

Fleet Numerical Weather Facility 
Naval Postgraduate School 
Monterey, California 93940 


Naval War College 
Newport, Rhode Island 02844 


Director, Naval Research Laboratory 
Attn: Tech. Services Info. Officer 
Washington, D. C. 20390 


AFCRL - Research Library 

L. G. Hanscom Field 

Attn: Nancy Davis/Stop 29 
Bedford, Massachusetts 01730 


Program Director for Meteorology 
National Science Foundation 
Washington D. C. 20550 


American Meteorological Society 
45 Beacon Street 
Boston, Massachusetts 02128 


Commander, Air Weather Service 
Military Airlift Command 

U. S. Air Force 

Scott Air Force Base, Illinois 62226 


Office of Naval Research 


Department of the Navy 
Washington, D. C. 20360 


169 


Copies 


No. Copies 


25. Department of Commerce, ESSA 2 
Weather Bureau 
Washington, D. C. 20235 


1/70 





UNCLASSIFIED 
Security Classification 


DOCUMENT CONTROL DATA - R&D 


(Security claseification of title, body of abetract and indexing annotation muet be entered when the overel! report ie claeeified) 


1. ORIGINATING ACTIVITY (Corporate author) Zea. REPORT SECURITY CLASSIFICATION 
UNCLASSIFIED 












Naval Postgraduate School 
Monterey, California 







3. REPORT TITLE 
NUMERICAL FORECASTING OF CLEAR AIR TURBULENCE 







4. DESCRIPTIVE NOTES (Type of report and inclusive detee) 


Thesi 
5. AUTNOR(S) (Leet name, fret name, initial) 


ETTEL, Michael J. 
MORGAN, William A. 


6. REPORT DATE Ja. FOTAL NO. OF PAGES 
. 170 
4 eli alii ; ; ! 


Se. CONTRACT OR GRANT NO. Se. ORIGINATOR’S REPORT NUMBER(S) 





b&b. PROJECT NO. 


Whis documd9?- ain ER be a NO(S) (Any other numbere that may be aeelgned 
*Welease -n| |; ” arc. £ 
! f Aves. or eae 


’ 


“ by ‘ + 


11. SUPPL EMENTARY NOTES | 12. SPONSORING MILITARY ACTIVITY 
Naval Postgraduate School 
Monterey, California 


13. ABSTRACT 
There is much disagreement as to (a) what causes clear air turbulence 


(turbulence which is not in or near convective clouds and is above 15,000 feet 
in altitude) and (b) which meteorological parameters can be used to detect and 
forecast its occurrence. The approach to this problem has been to relate not 
one parameter to clear air turbulence but various parameters. By summing these 
parameters areas can be defined where there is a high probability of encountering 
clear air turbulence. Each parameter has been based on a statistical study which 
found a relationship with clear air turbulence. The parameters used were hori- 
| zontal and vertical shear, curvature, kinetic energy and their derivatives. The 
numerical forecasting program proposed here can be extended to the stratosphere 
dwhen more reliable height and temperature fields are available. This program will 
+ have much more significance when intermediate forecast height fields, temperature 
i fields and a grid of much smaller mesh length are available. 





DD 59%. 1473 UNCLASSIFIED 


Security Classification 


171 





UNCLASS IF LED 


~ Security Classification 


14. (te 
REY WORDR LINK 8 
ROLE ROLE 


Turbulence 


Clear Air Turbulence 





wer geet eine 
> aaa’ S38 v= ' , a Aes pal, . 
= CN oat “ ey 8 iit a wl 
4h 7 ; ! 
1; 
are 
(BACK) 
1 
i NOV 68 UNCLASSIFIED 
S$ Oo © a 2 ee ; 
101-807-6821 Security Classification A-31409 























